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As a preliminary step, the radiation produced by a classical charged current coupled to 
a quantized is solved. To each order in a, all infrared divergences cancel between the 
virtual 7's and the real 7's absorbed from the plasma or emitted into the plasma. When 
all orders of perturbation theory are summed, the finite answer predicts a suppression of 
radiation with uj < aT. The analysis of QED then consists of two steps. First, a general 
probability at T 7^ is organized so that all the virtual e^, 7 are in the amplitudes and 
all the real e^,7 are in the phase space integrations. Next, the cancellations of IR 
divergences between virtual and real are demonstrated. 
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1. INTRODUCTION 

In studies ofthe quark-gluon plasma to be produced in ultrarelativistic heavy ion 
colisions and, more generally, in studies of field theory at finite temperature, a central 
concern is how the cancellation of infrared divergences comes about. Infrared diver- 
gences occur when an on-shell particle {p^ = m^) emits a massless particle {k^ = 0). 
The resulting propagator is [m^ — {p — ^)^]~"^ = [2p • k]~^ = [2{E — \p\ cos^)|^|]~-^ and 
so amplitudes for one- photon emission behave like ~ 1/k at low energy. At T = this 
behavior leads to logarithmic divergences, J dk/k, for each real and virtual photon. At 
T 7^ it leads to linear divergences, J dkuB/k, because of the Bose- Einstein function 
IT'S — l/[exp(/c/T) — 1]. This paper discusses how all the infrared 

divergences cancel even when T ^0. 

2. IR CANCELLATION FOR A CLASSICAL CURRENT 
2.1 Bremsstrahlung to Order a 

If a charged particle scatters while passing through a fixed-temperature plasma, it 
will radiate. If the radiated energy u is much smaller than the energy transfer in the 
collision, then the inelastic cross section factors: 

^ da ^ dP(q^) da 

so that the coUison cross section da/dq^ is independent of the photon energy- momentum 
(u3 — \k\). To first order in a the probability of radiating is 

dP 

2^^^ = E I^M ■ J'\' [1 + nB{uj)]/{2nf (2) 

pol 

When the radiated energy u is small, the current has a universal form regardless of the 
spin of the charged particle: 

V^_P^\ fc/2A 



Here A is a momentum cutoff that is necessary later. The radiation is mostly parallel 
to p or p When integrated over angles the result is 



^ = ^[l+n^(a;)]e-/^ A{p.p') = ^ 

dw U3 IT 



iln(i±^)-2' 



(4) 



where v is defined hy p ■ p' = m^(l — v^)"-*^/^. At large momentum transfer (Q ^ m), 
the behavior is ttA ?a 4a\ii{Q/'m). Except for the statistical factor ub, (4) is a classical 
formula. When to <^ T this predicts dPi/duj ~ AT/uj'^. This is totally unphysical 
because the energy radiated in low energy modes, below some £^max, would be infinte: 
duiuidPi/ du3 — oo. 
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2.2 Bremsstrahlung to All Orders in a 

To improve upon (4), couple the classical current (3) to the quantized radiation 
field Ay,. The generating function for all multi-photon amplitudes can be obtained 
by a Gaussian functional integration. From this one obtains the exact multi-photon 
amplitudes M.. The probability that any number n of real photons in the plasma will 
radiate a net energy a; is 

^ = f; I d$i...d$„5(A;° + ...A;0-a;)l J]|A1(A;i,...A;„)|2 (5) 

n=l ' pol 

^^i^^^^i^l)\^i^1)^^B{\%\)\ (6) 

This weights photon emission {k^ > 0) with the statistical factor 1 + ub, and photon 
absorption {k^ < 0) with the factor ub- 

Each amplitude M. is infrared divergent from closed loops of virtual photons. Each 
integration over the real photons is infrared divergent. The virtual and real contri- 
butions exponentiate to give 

/d^k 
-^-^J^{k)J^{k) ([1 + ubV^'' + UBe-^^^ - [1 + 2nB]) (8) 

The term proportional to l+ns represents stimulated emission; the term proportional 
to Ub represents absorption; the term proportional to 1 + 2nB represents the virtual 
photons (emitted and absorbed). At small k, ~ 1/k and ns ~ 1/k. Nevertheless, 
all infrared divergences cancel in (8) and give a finite result. It is possible to compute 
R{z) and to compute the Fourier transform (7). The final result is 

i;~'^2^'2^^' Att^t r{A) [it J ^ ' 

where A is the same function (4) as before. The most interesting feature of (9) is the 
appearance of two dimensionless scales: ^ <^ 1 and oj/T. If AttT <^ oj then 

dP A 

^«-[l + nB(a;)]e-'^/^ (^ttT « a;) (10) 
dw uj 

which agrees with the first-order result (4). However this does not apply at a; ^ T. At 
small energy, 

dP AT . 
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Naturally (10) and (11) agree in the region of overlap. Surprisingly, dP/dco is constant 
for u <^ AttT rather than increasing like l/u'^. The interpretation of the small u 
suppression is that the quantity 2AtvT = is a damping rate produced by the radiation 
reaction as required by unitarity [1]. 

3. REAL & VIRTUAL PARTICLES IN THERMAL FIELD THEORY 

We now set aside the semiclassical approximation and turn to the full quantum field 
theory. Each species of particle has four types of propagator in thermal field theory, e.g. 
Sab for electrons and D^^ for photons with a, 6 = 1 or 2. It is possible to rewrite thermal 
probabilities as squares of amplitudes that contain only and Df^^ , integrated over 
physical phase space. 

For definiteness, consider the process e~{pi) + plasma —>■ anything, with rate 

F,I 

Using completeness and thermofield dynamics [2] one can write this as 

i?(pi) =< m\c^c\m >= E I < > I' (13) 

F 

where |F(/3) > is a compete set of thermal modes. These correspond to real particles: 

oo „ oo „ ^ 

RiPi) = Yl / d^2...d^iYl / d^i...d^n-AMe,n{pi, ■■■P£;ki, ...kn)f (14) 

1=2 n=0 

where Aign is the amplitude for £ real e's or e's (including the initial one) and n real 
7's. By charge conservation, £ is even. Here is the fermion phase space similar to 
the photon phase space in (6). The explicit amplitude has the form 

M,,, = 5]</|[a#,[a#,...{6(p2),[-S,6t(p,)]}]]|/> (15) 

These amplitudes are related to Green functions that do not have the usual ex.p(—PH)/Z 
but rather S exp{—/3H)/Z. Consequently Mi^^ has no "cut" propagators D21 = D12 or 
= — 'S'i2- It is constructed entirely of Dii{k) and Sii{p). See also [3,4]. 

4. CANCELLATION OF INFRARED DIVERGENCES 

To analyze the rate (14) one can repeat the analysis of Yennie, Frautschi, and Suura 
[5] with slight modifications. Virtual photons inside the amplitude Aign can cause an IR 
divergence when they are on-shell and attached to "external" fermions pi,...pi. Each 
end of a photon line attached to an "external" fermion gives a multiplicative factor 
e2p^ I2p ■ k plus non-IR terms. At each order of a, the amplitude has a maximum IR 
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divergence plus non-leading divergences dP ^ ^ dP When summed to all orders, 

the leading and non-leading divergences all exponentiate: 

M,„ = exp(V^) MS ^ = / + ^^^^ e-'/A (10) 

V by itself is IR divergent. For real photons, infrared divergences occur in the integra- 
tions in (14). For fixed fermion momenta, define 

oo „ 

Rl{Vl.-Vl) - \ d^X-d^n-Mt^n? (17) 

n=0 ^ 

Since ~ kdk[6{k'^) + ub]-, divergences arise when Men ~ l//c- When n real pho- 
ton lines are attached to the external" fermions in all possible ways, they each give 
multiplicative factors e2p^ /2p ■ k plus non-IR terms. The contribution to (17) of n real 
photons has an IR divergence a" plus non-leading divergences q^~^ , a"^"^, .. 1. When 
summed over n, the leading and non-leading divergences all exponentiate to give: 

Ri{pi, ...pe) = Y du— d^i...d^rnP£,m{ki,...km)/m\ (18) 

All IR divergences are contained in the quantity 

This integration is IR finite because of the delicate cancellation between the virtual 
photons and the real photons (emitted -|- absorbed). The calculation is the same as in 
Sec. 2 with the result 

rIP A 1.) p'^/'^T -\u>\/K n rp 



but now A depends on all the "external" fermion momenta 

e 

87r2 Vab "A 1 - '"ab 

a, 0=1 

and Vab, defined by pa -pb/m? = (1 — 'Vab)~^^'^, is relative velocity of charge a in the rest 
frame of b. 



V ^-lnfi±^~) >0 (22) 
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5. CONCLUSIONS 

For fixed momenta of the "external" fermions, the rate (18) is IR finite because 
(a) the function dP/du is finite and (b) the integration J dwdP/dw is finite due to the 
behavior dP/dui — > constant as a; — > 0. 

The full rate requires integration over the thermalized fermions: 

CO „ 

^(^'i) = Y. d^2...d^i ReiPi, ...Pi) (23) 
e=2 

The fermion integrations do not affect the infrared finiteness of R{pi). However, there 
will still be Coulomb divergences that arise when any one of the momentum transfers 
vanishes: (po — Pb)^ — 0. It is known that the usual T = Coulomb divergence 
J sin 9d9/0'^ is reduced to logarithmic, J sin 9d6/ 9'^ at T ^ due to Braaten-Pisarski 
resummation [6]. 

For definiteness, the process e~ + plasma — > anything was treated specifically. 
For a general process {A} + plasma — ^ {B} + anything, where {A} and {B} are any 
sets of 6=*=, 7, the same arguments apply and show that all infrared divergences cancel. 
However, the logarithmic Coulomb divergence remains. 
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